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1 Abstract 

In previous galactic dynamo theories of the origin of the magnetic held in 
our galaxy, the subject of hux-freezing has been omitted. As a consequence, 
the equation of mass how has generally also been omitted, particularly in 
the halo where the galactic gravitational held will operate on the mass how. 
In this paper it has been shown that this neglect could have serious conse¬ 
quences for the results obtained from those galactic dynamo simulations that 
include the halo. A modihcation of these dynamo theories is proposed which 
involves the expulsion of very small bits of the magnetic held lines, rather 
than the wholesale expulsion of the complete magnetic lines encapsulated in 
the previous theories. This expulsion is accomplished by a spike instability 
that arises from superbubbles, when they break out of the galactic disc and 
their shell fragment. This leads to a cut in the lines of force that still remain 
in the disc. Subsequently, normal disc turbulence rotates the cut lines and 
thus dissipates their mean hux removing them from a role in the dynamo 
theory. This new process takes a length of time comparable, but a little 
longer than the previous growth time of the disc dynamo, but avoids the 
previous difficulties associated with flux-freezing and flux expulsion. 

2 Introduction 

In previous galactic dynamo simulations, the mean-field equations have been 
solved for the magnetic held, but the equations for the mass how have been 
neglected. This is not important for the solutions in the galactic disc, because 
there the mass hows are very slow. However, it is important in the halo in 
which the magnetic hux is dihused upward across the halo to entirely escape 
from the galaxy. This escape is necessary for our galaxy because the total 
hux must be conserved and the hux in any dynamo origin theory must start 
from a very small value and grow to its present value. In our galaxy the 
present hux is observed to be non zero. 

However, during the amplihcation the mass which is tied to the hux by 
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flux-freezing must also diffuse out of the galaxy with the magnetic flux. This 
mass is acted on by the gravitational held of the galaxy. This gravitational 
held can pull the mass and hux tied to it back into the galactic disc as fast as 
the halo turbulence is dihusing the hux upward. Indeed, a simple estimate 
shows that for the generally accepted value of the dihusion coefficient (of 
order of 10^® cm^/s ) the gravitational force dominates over the dihusion and 
the dynamo is strongly inhibited. This estimate comes from an elementary 
one-dimensional problem that includes dihusion, hux-freezing, and gravity. 

To get around this difficulty a new idea is presented in this paper for a 
process that dissipates the hux by the amount required by any dynamo origin 
theory. This process does not require the removal of mass from the galaxy, 
as is required by previous origin theories when hux-freezing is included. It 
involves those galactic superbubbles that break out of the disc. It is known 
that when a superbubble shell does break out of the galaxy its shells break 
up into many fragments, and from each of these fragments, a new instability, 
the spike instability can develop. This instability accelerates a small mass, 
and the small piece of the hux swept up by it, to faster than the galactic 
escape velocity allowing these pieces to ballistically escape from the galaxy. 
As a consequence, the part of the hux line attached to it remaining in the 
disc is ehectively cut into short pieces of hnite length. These pieces, although 
connected to each other through the expelled pieces, are not inhibited from 
rotating freely and in this way way flux of these lines is dissipated. Althongh 
very little mass is actually removed from the galaxy and the rest of the mass 
on the the line remains in the disc, its hux not longer connts in the balance 
of hux in the dynamo theory. It turns out that this in sit tu dissipation is 
equivalent to that resulting from hux expulsion, bnt it not snbject to the 
gravitational difficulty arising from the hux-freezing constraint. 

This paper is devoted to elaborating how this second process can enable 
the alpha-omega dynamo to fnnction without the difficulties involved in hux 
and mass escape. 

In the next section, an attempt is made to quantify the difficulties en¬ 
countered when hux-freezing and the ehect of gravity on the related mass are 
included in the dynamo simulations. These are semi quantihed by a simple 
one-dimensional model, and shown to be potentially serious. 

In the third section, I present the proposed new solution, and introdnce 
the spike instability and the relevant properties of the superbubbles. In sec¬ 
tions fonr and hve, I describe the snperbnbble phenomena and its properties 
which underlie the process. In section six, I derive the linear and non linear 
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properties of the spike instability. In section seven, I give an estimate for the 
number of spikes, the rate at which any relevant line is cut per unit time, and 
the mean length of the cuts. In section eight I calculate the rate of rotation 
of the relevant cut lines of force, which depend on their length, and the time 
for their flux to be dissipated by this rotation. This time is the important 
time which sets a limit on the speed of the dynamo if no line escape is pos¬ 
sible. The number of lines whose lengths suffer the spike instability is also 
estimated. 

The tentative conclusion from this: in approximately a billion years 
enough cuts are produced that every line in the outer quarters of the disc, 
loses its mean flux. This implies that the dynamo can grow the galactic field 
without the need for actual flux escape. However, the rate of dynamo growth 
is somewhat slower than thought at first. 

3 The inhibition of previous galactic dynamos 
by flux-freezing 

Our galaxy has a strong large-scale magnetic held and its origin is of consid¬ 
erable astrophysical interest. Where does it come from? It is generally agreed 
that the galaxy started off with a very weak seed held, perhaps as weak as 
10“^® gauss produced by cosmological processes, (Kulsrud et. al. 1997 ; Moss 
and Sokoloh 2013). This held is believed to be increased exponentially by a 
galactic dynamo. The problem is complicated for our galaxy, which is known 
to have a non zero hux frin measurements which indicate that the held does 
not change sign across the galactic midplane. This implies that some hux 
must leave the galaxy. 

A large number of important papers have been written on this galactic 
dynamo. Examples of these include the following references: Baryshnikova 
et. al. 1987; Brandenburg et al 1992; Brandenburg et. al. 1993; Dobler 
et. al. 1996; Donner and Brandenburg: 1990; Moss and Brandenburg 1992; 
Moss, et. al. 1993; Poezd et. al. 1977, Ruzmaikin et. al. 1988 Beck et. al. 
1996; and Parker 1997. In these papers the structure of the galaxy is taken 
as a thin disc with thickness 2H ^ 400 — 1000 pc. This disc is embedded in 
a large sphere of radius 4 — 10 kpc. The region of the galaxy outside of the 
disc is called the halo. 

In these papers the mean-held dynamo equations (Steenbeck et. al. 1966) 
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are successfully solved throughout the sphere with solutions that grow expo¬ 
nentially in a growth time of a few hundred million years. 

These equations are 

dBr d{aBe) d Z' dBA 

dt ~ dz ^ dzy dz J ^ ’ 

dBe dQ d (^dBe\ 

-w = + 

In these equations alpha, which produces the famous alpha effect, turns 
toroidal flux into poloidal flux. f2(r) is the differential rotation of the galaxy, 
and it turns poloidal flux into toroidal flux. (3 is the coefficient of turbulent 
magnetic diffusion which mixes up the magnetic helds. To allow growth of 
the galactic held, which starts with zero hux early on and at present has non 
zero hux, conditions must be satished on the surface of the halo that allow 
the escape of hux into the intergalactic medium. 

These integrations are mathematically correct, but they ignore two im¬ 
portant physical ehects, which throws some doubt on the physical soundness 
of their results. 

These two ehects are: 

1. Flux-freezing in detail. The mean held equations given in Equation ([1]) 
do preserve hux-freezing on the large scale, but they do not take into account 
that even on small scales any mass tied to a given hux tube must remain tied 
to this hux tube as it moves. By ignoring this property the solutions have 
been able to avoid the dynamics of the mass how and are solved without 
regard to the mass dynamics, which would follow from hux-freezing. 

2. The gravitational held. Once mass how is ignored there is no need 
to take into account the gravitational held of the galaxy. However, it has 
a strong ehect on any mass how tied to the hux how. In any simulation in 
which the mass is tied to the hux, the gravitational held plays an important 
role. 

In order to understand the importance of these ehects most clearly, I 
consider only the earliest phases of the galactic dynamo. This is the phase in 
which when the Lorentz forces are too weak to be macroscopically important, 
but strong enough to enforce hux-freezing. 

I show by a simplihed model, that when these two ehects are included in 
the simulation, the gravitational held will strongly ahect the how of the hux 
out of the halo. To this purpose, I ignore all dynamo terms except the beta 
terms, and I assume that (3 is constant. 
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The flow of the flux is treated most simply by considering what happens 
to a bundle of flux which leaves the disc at time t = 0. I assume one 
space dimension, z, in the vertical direction normal to the disc. Let / be 
the intensity of the bundle and let us describe it by a collection of particles 
diffusing with diffusion coefficient f3, and under the influence of a downward 
gravitational field g. It is physically reasonable to assume that the diffusion 
will occur in an downward accelerating frame with a diffusion coefficient 
whose value is independent of the acceleration. The frame is accelerating 
downward with at a constant value g. 

Then as shown in the appendix, $(t) the fraction of this bundle of parti¬ 
cles that is beyond the radius of the galaxy, R at time t is. 


$(t) = Erfc \{R + gf/2)/^/4^t 


( 2 ) 


where Erfc is the compliment of the error function. 

If 5 ^ = 0 the simplified model corresponds to the previous simulations and 
the fraction escaping is simply 


m 


Erfc 


R/VWt 


( 3 ) 


As t increases, the argument of the complimentary error functions de¬ 
creases, and <h increases. Finally, when when tc where 

r2 

R = — = 8 X 10^i?io7/928 years (4) 

the argument of Erfc is much less than unity , and the bulk of the flux 
is beyond R. Riq is the radius of the halo in units of 10 kpc and /328 is jS 
in units of 10^^ cm^/s. These values are consistent with the dynamo times 
times attained in the dynamo models, which shows that the example gives 
roughly the correct results. 

If we include the galactic gravitational field, g ~ 10“® cm^/s, then the 
fraction of lost particles (or flux) is given by the full Equation ([2]). In this 
case, as t increases, the argument of Erfc at first decreases till it reaches a 
minimum value at t = tg, where 


t 


9 ~ 



( 5 ) 
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and then increases. 

That is, the lost flux in our model increases to and then seems 

to decrease. In our model this corresponds to particles that have crossed R 
and then fallen back inside z = R. Presumably, conditions beyond z = R 
are different. For example, g may weaker and other factors may enter not 
included in the model. Therefore, assuming any flux that has passed beyond 
the halo boundary escapes, the largest fraction of the flux that can escape is 




Erfc[ 


(2i?/3)3/4^V4 

^1/2 J 


( 6 ) 


tg is essentially the free fall time from the radius R. Numerically, we have 
the following table. 

For i? = 10 kpc and (3 = 10^®cm/s, <F = 

For R = 5 kpc and (3 = lO^^cm/s, <F = 

For i? = 10 kpc and (3 = lO^^cm^s, <F = 0.016. 

For R = 5 kpc and {3 = lO^^cm^s, <h = 0.191. 

Thus, we see that if [3 is of order 10^®cm^/s (Poezd 1993, Dickey and 
Lockman 1990, Kulkarni and Fich 1985) the fraction of the flux that escapes 
is extremely small, and that nearly all of the flux escaping from the disc is 
reflected back by the halo. That is to say the disc boundary conditions are 
essentially closed conditions, significantly modifying the standard dynamo 
solution. 

This simple approximate model implies that at the very least the standard 
simulations should be repeated with mass attached to the flux, by flux¬ 
freezing, and a gravitational field included. This should be the case even 
when one has a larger diffusion coefficient. 

by the way, if the diffusion coefficient were large enough, (i. e. 10^® cm^/s) 
to produce the escape of the flux in the presence of flux-freezing and gravity, 
a considerable amount of energy is involved. In fact, if a sufficient amount of 
flux escapes to allow the dynamo to double the magnetic flux in the disc, at 
least half of the disc will be lifted out of the galaxy. This will take approxi¬ 
mately 10^^ ergs for an annulus with width a kpc. If the removal occurs over 
a time of a billion years, the amount of power that must be supplied is of 
order 10^® ergs per second. This is and interesting number. The power need 
to resupply the loss of cosmic rays is a little larger. About t times 5 x 10*^° 
ergs per second and infalling mass has been invoked to supply this energy. 

Clearly, the above argument rests heavily on the property of flux-freezing. 
To see that flux-freezing is valid in this case, consider the magnetic diffusion 
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( 7 ) 


equation, in a static medium, with non turbulent resistivity. 


^ = ilWB 

dt in 


where according to Spitzer (1962) rj/in = 10^/Te'^^ civ?/s where Tg is the 
electron temperature in electron volts. Let B have a scale length of ^ and 
estimate ~ B Then with Tg = leV (a typical astrophysical temper¬ 
ature) , 


5B Wt 3 X 1023 

- - 

B P e 


3 X 10"^^ 


( 8 ) 


if t = 10® y and ^ ~ 1 pc. If the plasma is moving, then this change in 
B represents the slippage of magnetic flux relative to the plasma. For the 
galactic dynamo these are appropriate values for t and ^ as the flux leaves 
the disc, and this very small slippage justifies the flux-freezing assumption. 
In any event, if [3 is not sufficiently large, then either the standard galactic 
dynamo and needs to be revisited or flux-freezing must be invalid. 

In light of this, it would seem prudent to look for a different way to get rid 
of the flux, and it is the purpose of this paper to suggest such an alternative. 


4 Alternative suggestions for a solution to 
the problem 

Clearly, if the halo diffusion coefficient is not large enough something more is 
needed to make the galactic dynamo viable. Parker (see hgure 3c in Parker 
1992) has suggested that it is not necessary for entire lines to be bodily moved 
out of the halo, but only partial loops of the lines. The mass on the loops 
themselves could be greatly reduced by mass sliding down the held lines, 
and the loops lifted upward by buoyancy or cosmic ray pressure. Parker 
then suggests that these loops would be cut loose by magnetic reconnection 
and lifted out of the galaxy by buoyancy and cosmic ray forces. However, 
Parker’s lifting mechanism is not available in the early very weak held phase 
of the galactic disc. It is also true that in Parker’s model most of each hux 
line would be left behind and so the disc hux would appear to be unchanged. 
However, in section 9, I show that, even without reconnection, such pieces 
of the held lines could be dissipated by horizontal turbulent dihusion. 

Other plausible approaches to the hux removal problem are the use of 
supernovae or superbubbles to lift portions of the held lines out of the galaxy. 
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In these approaches, the main flux is also left behind in the disc, although 
these remnants of the held lines might also be dissipated as in section 9. 
However, neither cosmic rays, supernovae or superbubbles have the capacity 
to directly lift even pieces of enough lines out of the galaxy. Most type 
II supernovae are too embedded in the disc to plow their way through the 
interstellar medium and leave the disc. Only the very few type II supernova 
that explode high above the disc could throw material out of the halo, and 
there are not enough of these. The largest superbubbles potentially have 
enough power, but once their shells leave the disc they become unstable and 
break up into fragments which fall back into the disc (Raflkov and Kulsrud 
2000, Mac Low et. al. 1989). 

This is true of our galaxy, but not necessarily true of others. In fact, 
in starburst galaxies the sources of activities are much stronger and strong 
plasma outflows from then are observed. But our galaxy is not forming stars 
at the rates of these galaxies, and there is no evidence that it ever did. On 
the other hand, smaller galaxies have weaker gravitational fields that could 
allow outflows. 

Thus, the possible solutions to the flux problem in our early galaxy face 
two difficulties: (1) First, the existence of enough sources sufficiently powerful 
to expel mass from the halo. (2) Second, the removal of the net flux of any 
piece of a line left behind. 

In this paper I suggest a different approach to the problem: 

First, I propose that fragments of a superbubble shell which breaks out 
of the disc are unstable to forming spike-like structures that can drive small 
lumps of mass upward. The mass in these spikes acts like the mass in Parker’s 
loops, in that, mass can slide down their sides reducing the mass on their 
tops. The mass remaining on their tops can then be accelerated to the 
galactic escape velocity, and escape ballistically out of the galaxy. This mass 
will carry a small portion of the lines of force out of the halo. This will 
resolve the first difficulty. 

Second, it is true that the line removal leaves the rest of the line in the 
disc, but this remainder is effectively cut into separate pieces. These pieces 
will be rotated by the turbulent motions in the disc. These pieces are still 
connected to each other through the pieces of flux expelled by the spikes. 
However, it easily seen that the flux in the disc is thus reduced, the residual 
part of the flux being carried away by the tiny pieces of expelled flux. Thus, 
surprisingly, the flux will be reduced by the cuts due to the spikes, without 
the bulk of the flux lines or the mass being removed from the disc. 




Figure 1: The line breaking with and without the spike 


If there were no spikes the lines of force could abstractly still be considered 
to be divided into pieces much as they are with the spikes. Because the 
magnetic held is very weak the pieces could still rotate randomly. But in this 
case, these pieces are connected to each other through hux lines still in the 
disc so that no hux is lost. 

The two cases are illustrated in hgure 1. On the left we have a line of 
force AD in the xy plane. Imagine it divided into two parts AB and CD. 
Let the two parts be rotated, as in the hgure, to A’B’ and C’D’. These two 
parts have no hux in the disc in either the x direction or the y direction . 
However they were connected by a short line BC which by hux-freezing goes 
to B’C’, so that the combined line A’B’C’D’ has the same hux in the disc as 
the original line. 

Now, consider the same situation and the same line segments, except B’C’ 
has been replaced by a line that stretches out of the disc to the spike and 
back. This segment has very little hux in the disc and the remaining lines 
have no hux. As far as the disc is concerned the line has been cut into two 
parts whose disc huxes can cancel. 

It should be noted that if the line B'C reconnected with itself, then the 
matter along it would naturally drop back into the disc and draw the line 
B'C with it. Thus, reconnection would destroy the model. But it seems 
highly unlikely that hux-freezing on this scale would allow reconnection. 
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As the cut pieces rotate and decrease the disc flux, the magnetic line in 
the expelled parts will twist in the opposite direction consistent with flux 
and helicity is conservation. However, the helicity in the disc, under any 
plausible deflnition of local helicity, will change, as required for any dynamo 
that is growing the held. 

In the rest of this paper I show that in a billion years, there are enough 
spikes produced by superbubbles to cut enough lines to reduce the flux by 
a finite value in the top and bottom parts of the disc. This will resolve the 
second difficulty. 


5 SuperBubbles 

Next, we discuss the superbubbles. 

Superbubbles were first discovered by Heiles in 1974. (See Heiles 1974) He 
found them by surveying the galactic disc in narrow, 10 km/s, width Doppler- 
shifted velocity bands of the 21 centimeter line. He found that in each band, 
the intensities were located in definite stripes. Combining these stripes in 
velocity and space, he interpreted them as gigantic expanding bubbles. These 
bubbles were similar to those produced by stellar winds, (Weaver et al 1977), 
but on a much greater scale. 

It was soon realized that superbubbles must be produced by multiple 
sequential supernovae originating in a forming galactic cluster of stars. The O 
and B stars in the cluster explode as supernovae almost as soon as they form, 
so that their energy instead of going into the kinetic energy of their remnants, 
goes into a bubble of high-pressure low-density plasma. This plasma pushes 
all the surrounding mass into a massive shell. An example of the larger of 
these superbubbles is the North Polar Spur. It has an enormous diameter of 
more than 80 degrees stretching across the whole sky. After a period of time, 
it was appreciated that nearly every superbubble could be identified with a 
young galactic cluster, and conversely every galactic cluster was generating 
a superbubble, (Elmegreen and Clemens 1985). 
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In these superbubbles the supernovae explode as frequently as once or 
twice every million years. The rate of explosions multiplied by the supernova 
energy per explosion gives the luminosity, L, of the superbubble, which ranges 
from 10^® to 10^® ergs per second. A superbubble lasts as long as the O and 
B stars in its cluster are forming. Generally, this lifetime is 25 to 50 million 
years. The average number of supernova in a cluster lies between 6 to several 
thousand. The luminosity of the superbubble L is emitted continuously and 
can be estimated a.s L SEgn, where S is the supernova rate and Egn is the 
energy of a supernova. 

It is found that the shells of superbubbles with L > 10^^ ergs/s, expand 
long enough to break out of the disc. Prior to breakout, these superbubble 
decelerate, as more and more mass is snow-plowed into their shells. However, 
for the stronger ones after they breakout of the disc, there is no further sur¬ 
rounding mass to snow-plow, and the pressure of the bubble then accelerates 
the shell. This situation is unstable to the Rayleigh-Taylor instability because 
the low density plasma of the superbubble accelerates the high density shell. 
As a consequence, the part of the superbubble extending out of the disc, 
fragments into many pieces. This is observed in the numerical simulations 
of Mac Low et. al.(1989). After the breakup, the fragments continue to rise, 
but eventually their mass is brought to a stop by the galactic gravitational 
held, at heights of one or two kiloparsecs. 

This does not directly help the dynamo because when the superbubble 
mass falls back into the disc it drags the hux back, thus expelling no net 
hux. However, it is shown in section 7 that a spike instability develops on 
any fragments that manage to hold onto the superbubble pressure for a long 
enough time. This spike will manage to throw mass out of the halo even if 
the fragment falls back into the disc. 

In fact, the mass at the very top of the spike achieves velocities greater 
than the galactic escape velocity. The mass at the top of the spike will thus 
drag the flux lines that are still embedded in it, out of the galaxy. Each 
of the very short segments of these lines will ’cut’ the rest of the line, that 
remains in the disc. 

The part of the shell above the midplane of the disc will contain all the 
flux in the top half of the disc and within a sphere of radius H. The spikes 
are unstable only on the top of half of the fragments and contain flux and 
mass from one quarter of the disc. Similar things happen to the part of the 
shell below the midplane. 

The rest of the paper, is devoted to Ending the number of cuts in a given 
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period of time. We first estimate the number of fragments into which the 
superbubble shell breaks. We assume that / of these fragments will have a 
spike. Then, after investigating the nature of the spikes, we calculate the 
number of lines of force it ejects from the galaxy, i.e. the number of cuts it 
makes. 

Next, we hnd the total number of spikes and cuts above the midplane 
due to a single superbubble Then, making use of the rate of production of 
superbubbles per unit area and time, we calculate the number of cuts made 
per square kiloparsec per million years. 

By considering an annulus through the solar circle, we determine both 
the number of lines and the number of cuts in this annulus in fg million years. 
The ratio of these two numbers is the number of cuts per line in this time. 

From this number we get the average length of the cut lines at any time. 
From diffusion theory with the a diffusion coefficient f3d obtained from dy¬ 
namo theory, we calculate the rate of spreading of 6d the angle of the cut 
lines relative to its initial direction. The time it takes for this average angle 
to become larger than 7r/2 is the dissipation time of the flux. On comparing 
with the standard value for the dynamo growth time we can test whether the 
dissipation time is shorter or longer than the dynamo time of the models. 

The theoretical and observational properties of superbubble can be found 
in the following references. (Heiles, 1974, 1979, 1984, 1987, 1990, Weaver 
et al.l977, Mac Low and McCray, 1988, Mac Low et al 1989, McCray and 
Kafatos 1987, Elmegreen and Clemens 1985, Kennicut et al 1989, and Rahkov 
and Kulsrud 2000). 

6 Superbubble Properties 

For simplicity, we assume that the galactic disc is a slab of uniform density, 
no = 1.0 cm“^ and has a half thickness oi H = 200pc, (Ferriere, 1998). This 
slab contains the same interstellar mass as the galaxy. We assume that all 
superbubble explosions are centered on the disc midplane. From now on , 
we consider the parts of the shell and its fragments only in the upper half of 
the disc. 

For this case, the radius, R, of a single superbubble before it breaks out 
of the disc, evolves as 

R = 267Ljf 4 /^pc 

where L, the luminosity of the superbubble is lO^^Lsg ergs/s and t = t7l0^years. 
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(Rafikov and Kulsrud, 2000 ; Mac Low and McCray, 1988; Weaver et al. 
1977.) 

Differentiating this with respect to t we hnd the velocity of expansion 

V = 15.7L38^t7 ^'^^km/s 

From these results we can determine the luminosity Li at which a su¬ 
perbubble is just able to break out of the disc. We use the fact that the 
superbubble stops expanding horizontally when its velocity drops to the in¬ 
terstellar velocity of 10 km/s, and set its horizontal radius equal to H at this 
time. (Ferriere 1993) (Under the uniform density approximation for the disc, 
the superbubble will expand spherically up to breakout. Afterward, it will 
expand only in the vertical direction.) 

We hnd that 

Li = 1.46 X 10^'^ ergs/s 

and the breakout time for a superbubble with L > Li is 

tbr = 11.7 f — j Myr 

where Myr denotes 10® years. 

The rate of production of superbubbles is 

da = 0.10 ^kpc-2 Myr-^ (9) 

\Li J Li 

where we have normalized the rate to the Li luminosity. (Ferriere 1993, 
Kennicut et. al 1989, Elmegreen et al. 1985.) 

We are only interested in superbubbles with luminosities great enough 
to break out of the disc, i.e. L > Li. The shells of these superbubbles will 
fragment by the Rayleigh-Taylor instability, when R > H and they leave the 
disc. The growth rate of the instabilities with horizontal wave number k has 
a maximum of 7 = ^/kgj2 at kD = 1, where D is the thickness of the shell. 
{g ~ 10 “® is the vertical gravitational held produced by its stars in the disc.) 
(Ferriere 1998) 

We can obtain the thickness of the shell at breakout by comparing the 
pressure p in the hot core of the superbubble with the pressure UsKbT in 
the shell. The core pressure can be obtained from the energy deposited by 
the exploding supernovae into the core plasma during the breakout time. 
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{b/ll)Ltiyr. (The factor of 5/11 takes into account the work done by the core 
pressure accelerating the shell.) Dividing this energy by the spherical volume 
of radius H, gives 3/2 times the core pressure. The result of the calculation 

is 

n, = 4O.O(£2/^0-^)cm-=^ 

where 



The temperature is not known, but it is believed to lie between 100 and 
1000 degrees, so we normalize it to 300 K and set T = 3006*. The shell sweeps 
up all the interstellar matter in the sphere above the midplane, so 

D = = 1.66(6-2/30) pc ( 11 ) 

3ns 

The superbubble shell will break into fragments with radii on a scale of 
order Aq = 27r/A; = 27rD, with k = 1/D the wave number at maximum 
growth rate, so 

Ao = 10.4(6-2/30) pc 

Mac Low et. al. (1989) have simulated the explosion of a superbubble 
with L ~ L 38 , and have displayed the fragmentation, at the time of breakout. 
For this luminosity, Aq = 6 pc. A crude inspection of their hgure 2c shows 
that the radii of their fragments at breakout is between 12 and 15 pc. We 
use this to set the radii of the fragments to 2 ^Ao, where ^ is of order one. 

The number of fragments per superbubble is thus 

Af/ = = 91.7(C'“ry-'") (12) 

Each fragment will have a spike with probability /. 


7 The Spike Instability 

Assume a thin disc at the top of a fragment of the superbubble shell with 
surface density (Jq. Ignoring the inertial forces, the gravitational force on a 
unit area of the shell is agg. g = IQ-® cm s~‘^, the real galactic gravitational 
held. 
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Figure 2: The area for the spike Instability 

Assume that the spike perturbation is azimuthally symmetric about some 
vertical axis. We derive the Lagrangian equations for it, with radial and verti¬ 
cal Lagrangian displacements R{rQ, t) and Z^tq, t), where tq is the Lagrangian 
radial independent variable. 

In the perturbation, consider a thin zonal disc generated by the rotation 
of a line initially extending from tq to tq + S. After a time t these two 
end points are displaced to [rg -|- R{rQ,t), Z{rQ,t)] and [rg -|- 5 -|- R{rQ,t) + 
SR'{ro,t), Z{ro,t) + (5Z'(ro, t)], where the prime denotes a derivative with 
respect to tq. 


If we rotate this displaced line by a small angle e about the cylindrical 
axis, we will get a small vectorial area A whose coordinates are 

Ar = -eS(ro + R)Z' (13) 

Az = eSi^TQ -\- R){i -\- R') 


See hgure 2. 

Now m, the mass of this area, is constant under the Lagrangian displace¬ 
ment, and equals derocro. The pressure force has the components pAr and 
pAz, with p = aoQ also a constant. The gravitational force has the compo¬ 
nents [0,—mg]. 
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Therefore, the (nonlinear) Lagrangian equations of motions are 


= -^{ro + R)Z' (14) 

ro 

= —i'i^o + R)i^ + R')~9 
ro 

Take R and Z proportional to e'*'*. Then the linear equations are 

j^R = -gZ' (15) 

7"z = <A- + H) 

ro 

Setting 7 ^ = kg^ where k gives the initial spatial dependence of the mode, 
we hnd the linear solution is 


d^R 

dR 

d^Z 

dR 


Z = ae^*Jo{kro) (16) 

R = ae^^Ji{kro) 


where Jq, andJi are Bessel functions. 

When the hnite thickness of the shell D is taken into account the growth 
rate is given by 7 ^ = kg /(1 + k'^D'^) so that the maximum growth rate occurs 
when kD = 1 , and is 7 = ^Jg/2D. Using equation (ITT]) for D we hnd the 
linear growth time of the spike is 7 “^ = million years. 

We need the nonlinear development of the solution to study the escape 
of mass out of the galaxy. 

The nonlinear limit of our equations is 


d^R 

dR 

d^z 

dR 


To solve these equations, set 


R 

Z 


-—RZ' 

ro 


-RR' 


ro 


giro) 

{to - ty 
i^{ro) 
{to - ty 


(17) 


(18) 
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Then the nonlinear equations become 


Qrj = (19) 

R 9 / 

bu = —rjri 

ro 

Canceling t] in the hrst equation, we get 

u = C-3^ 

9 

where C* is a constant to be determined. The second equation gives 

2 eC 2 9 4 12^ 2 

The linear and nonlinear solutions do not overlap. To a sufficient approxi¬ 
mation, we can compare the two solutions for Z at t = 0, to get rough values 
for C and to. 


C 3rg/ g 2^2 


From this we hnd that 


''0 


r - 

^r\ — 


c = 


a(l — k^r^/A) 


12 

k'^ga 


atl = 


( 20 ) 


( 21 ) 


12 12D^ 


k‘^g g 

Since kD = 1, the initial radius of the spike is \/2D and 

C = 3.14 X lO^®(r^/^0^) cm^s (22) 


Note that C does not involve the initial amplitude a. 

As t approaches to, both Z and R diverge. Before this happens, the 
velocity reaches the escape velocity and the top of the spike moves ballistically 
out of the galaxy. Setting the upward velocity dZ/dt at ro = 0 and t = tesc, 
equal to v^sc = 400 km/s, (Binney and Tremaine 1987) we have 

dZ _ 24 

dt k‘^g{Aty 
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so taking k = 1/D we have 


At = to — tesc = 0.37(£ Myr (23) 

The height, hesc, at which this happens is 

Aesc = = 75(r'‘''»»2/3) p, (24) 

At the same time the plasma at the top has spread out, and thus, the 
magnetic held originally trapped on the top will also have spread out. Its 
strength is weakened by a power of 


R 


IQC 1 

TW 


(25) 


at the time tesc- 


8 The Rate of Cuts 

As a result of a single spike, a line has been cut. See hgure 3. 

We will hrst estimate the number of cuts made by a single spike, then 
the number of cuts per superbubble, and then the rate of cuts in a cubic 
kiloparsec. 

How many cuts does a single spike make? From the previous section we 
see that the top of the spike will be accelerated to the escape velocity, and 
leave the galaxy, dragging its embedded hux with it. These lines will still be 
connected to the disc and each removed line will produce a cut in a disc line. 

Let us take the initial horizontal diameter of the spike to be d = 2D. 

Now, the lines in a spike will be those that come from the top half of its 
fragment, which in turn come from the top half of that part of the initial 
galactic held above the galactic midplane. In general, we do not know the 
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Halo 


After fall back 



With no tension the negative field diffuses away 


Figure 3: The Model for Line Cutting 
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vertical profile of the structure of the disc held (which is determined by the 
dynamo theory). But it must have the same prohle as the initial hux in 
the spike. Let the initial hux in the spike be xBo{‘2,DH/2) where the prohle 
factor X would be one if the held were uniform. 

Thus, the initial number of lines in the spike is 

0 ' = xBoDH (26) 

As the spike grows, the mass on the top slides downward along its sides, 
and at the escape time, the surface density is reduced by a factor of where 
p is given by Equation fl2SD . 

Although the surface density decreases by the number of held lines 
decreases only by one power of p. This is because in one direction, the 
mass hows along the held lines. The lines are only removed from the top by 
motions in a transverse direction to the held. Remember we are dealing with 
extremely weak helds, so the velocities are axisymmetric and the helds have 
no ahect on them. 

From Equation fl25|l we see that in the planer single-spike solution, the 
spike will expand by approximately thirty times. It is unlikely that the 
solution is valid to such a large distance, but the thinning out of the hux 
and density about the origin is locally correct. Since I am uncertain what 
terminates the nonlinear solution in the horizontal direction, I conservatively 
assume that the hux is ejected from the galaxy at least a distance tD from 
the origin, where r is a factor greater than one. 

Thus, the hux which is removed from the galaxy by a single spike and 
the number of resulting cuts is 

, tDHB DH^B ^ 9 . ^ 

‘^spike = ^—~— = ^ gauss cm^ (27) 

Next, the number of cuts per superbubble is the number of cuts per spike 
times the number of fragments times the probability / that each fragment 
will have a spike. Thus, 

= fNfct> gpike = 0.0236 x xB‘^B{tl'^B-‘^l^ (28) 

Now, using (Equation ([9]) the number of cuts per unit area summed over 
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all L, is 


0kpc J 

/ OO 

= .OO276xi/^5(0-2/3/r^-2) kpc'^ Myr”^ 

To get the number of cuts per unit time of any individual line in the 
upper quarter of the disc, we consider an annulus of radius 8.5 kpc, thickness 
H and width /S.R. The rate of cuts in this annulus is 

0ann = 27ri?s'Ai?0kpc (30) 

= 2.76 X lO“^xi7^5(27ri?5Ai?)(0-2/V/r^) Myr"^ 

The initial number of lines in the top quarter of the disc and in this 
annulus is 


T = 0.5xBHAR gauss cm^ 

Clearly, the total of cuts which any line in the annulus suffers in 

R Myr is 


^cuts = ^ = 1.10 X 10-^{27rRsMe-^/^Tfr^) 

The prohle factor, y has dropped out of this ratio. 

The average length of every line segment after tg Myr is 


27tRs 

^cuts 

= kp. 


(31) 


9 The Dissipation Time for the Flux 

We now consider how fast the cut lines can rotate in order to remove the 
net flux due to them. This rotation occurs in the disc and is produced by 
turbulent diffusion there. The disc diffusion coefficient is included in all 
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previous work. Is has usually been estimated from the observed properties of 
interstellar turbulence, but the magnitude of this turbulence is uncertain. In 
the meantime Ferriere has carried out a more deductive derivation of it that 
is produced by the fluid motions produced by supernova and superbubble 
explosions. (Ferriere 1993). Characterizing the effective turbulence by the 
disc diffusion coefficient 13d- she hnds v f3d ^ 5 x lO^^cm^/s. This diffusion 
coefficient is much smaller than the halo coefficient, discussed in section 3, 
and is the same as the beta in the dynamo equations (|Ti) . 

Consider any given cut portion of a line aligned along any axis. Each 
end will be displaced by the interstellar turbulence, a perpendicular amount 
whose standard deviation is 2f3t. Therefore, the line will be rotated from an 
initial direction by an angle 6d whose standard deviation is 

< {0dY >= l^ = (32) 

where tg is the time in 10 ® years. 

The dissipation time, td, is found by setting < 6 *^ >= ( 7 r/ 2 )®. This is the 
time when all the relevant flux is diffused to zero. 

td = 1.66 X lo9(04/9y-2/3^-2/3^4/3) (- 33 ^ 

This estimate for td shows that the time to dissipate the flux is of order 
of a billion years. In deriving it, I make use of certain assumptions that are 
quantihed by the parameters: 

• /, the probability that a fragment will have a spike instability, which 
is equal to the probability that the fragments holds their superbubble 
pressure for a time tesc, [see Equation fl23ll ]. 

• the size of the fragments normalized to AnD. 

• 6, the shell temperature normalized to 300 degrees Kelvin. 

• r, the radius of the top of the spike at the escape time normalized to 
2D. 

These parameters are included in Equation (|33ll . Their values are uncer¬ 
tain and and represent its uncertainty. In particular, if the spike instability 
did not occur at all (because all the fragments lose their superbubble pressure 
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too rapidly), then / would be zero, Equation for would be infinite, 
and the spike hypothesis of this paper would not work. 

For comparison let us take an estimate of the dynamo time tdynamo equal 
to the growth time, from one of the referenced papers in section 1, (Bran¬ 
denburg et. ah 1992). 

tdy = ~ 3 X 10®years (34) 

a time considerably shorter than a billion years. Thus, if the Brandenburg et. 
ah galactic dynamo is inhibited by flux-freezing and the galactic gravitational 
held, it still can work using this new process, but with a longer growth time 
given by td- 

Although we cannot determine the values of the parameters, we can make 
certain remarks. First / is the most uncertain parameter since the analytic 
details of the fragmentation are not known. / must be less than one, but 
is probably bigger than one third. The parameter r is dehnitely larger than 
one but less than since is the radius of the fragment. It is probably 

considerably larger than one. One might assume that the factors / and r 
cancel. ^ which has been determined by comparison with the actual simu¬ 
lation of Mac Low et. ah (1989) is probably not far from one. Thus, it is 
reasonable to believe that td is directly proportional to 6. 9 can easily be 
as small as 1/3 if the shell temperature is lOOiL or as large as three if the 
temperature lOOOK so this probably gives a plausible range for td- Thus, one 
might have either td is comparable with tdy in which case the previous values 
for the dynamo time are believable, or it may be much larger and the dynamo 
would be correspondingly slower. Some detailed numerical simulations are 
probably necessary to decide which is the case. 

How does one analytically combine our result for td with the alpha-omega 
differential Equation ([1])? The simplest way is to modify these equations by 
adding terms —Br/td and —By/rd to the right hand sides of the two equa¬ 
tions when .5H < \z\ < H. In addition, one must include the mass flow 
expressed in terms of the flux flow from the flux-freezing condition, and the 
gravitational held. It should be acceptable to use vacuum boundary condi¬ 
tions on the edge of the halo. This leads to a modihcation of the standard 
approach and should be more realistic. 

The broken pieces of hux from one growth phase will not interfere with 
the next growth phase because the strength of their helds are relatively small. 
This only holds during amplihcation. When the stronger held is approached. 
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the spike instability saturates and one can no longer neglect the residual 
held in the disc even after it is dissipated. This will probably saturate the 
dynamo as well, since td would be much longer, and the boundaries would 
become completely closed. Assuming saturation at the present epoch these 
loose finite pieces could be responsible for the random helds measured by 
Rand and Kulkarni (1988). 

Incidentally, we can calculate the length of the cut lines averaged over a 
decay time. This is given by averaging Equation (1^ over td, and is equal 
to the length at td, multiplied by a logarithmic factor which is two or three. 
This length of only a few kpc is about the bulk distance cosmic rays can 
propagate during their lifetime. These cut lines remaining after saturation 
suggest a natural method for cosmic rays to escape from the disc. 


10 Summary and Conclusions 

I have shown by a simple example that, when flux-freezing and the gravi¬ 
tational held are present the rate of expulsion of hux in previous galactic 
dynamo theories, can be strongly inhibited by the mass on the the outhow- 
ing hux falling back from the halo into the disc. A crude estimate of this 
inhibition is given by Equation (|2]) which indicates that if the halo dihu- 
sion coefficient is much less than 10^® cm^/s the outhow will be completely 
inhibited. 

However, in this case the dynamo may be saved by another process based 
on a spike instability arising from the fragmentation of those superbubble 
that succeed in breaking out of the disc. This instability can expel a very 
short length of a line of force, cutting the rest of the line in two. If a line is cut 
into short enough pieces they can rotate horizontally in the disc, dissipating 
their net flux and removing it from playing any role in the dynamo. This 
turns out to be equivalent to expelling the line as far as the dynamo is 
concerned. An estimate of how many cuts occur shows that the relevant 
amount of flux can be removed in a time of order of a billion years. Thus, the 
dynamo can still operate, but on this longer time scale, to grow the galactic 
magnetic held. The two processes, the alpha- omega dynamo and the hux 
annihilation by cutting helds, can be combined in a systematic theory. One 
merely has to add terms giving the annihilation rate, to dynamo equations, 
and assume boundary condition the compatible with the inhibition of the 
hux expulsion. 
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Incidentally, the method of dissipating flnx withont removing mass may 
apply to other dynamos, snch as the solar dynamo, in which the flux is 
conhned by flux-freezing, and one does not want to remove an appreciable 
fraction of the mass of an object to produce its dynamo. 

Although no direct observational evidence is produced for this model it 
may have indirect consequences of importance such as an explanation for why 
the observed random magnetic helds are larger than the mean held. Also the 
cut helds may provide an escape route for cosmic rays. 
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Appendix 


13 The Diffusion of Flux in the Presence of 
a Gravitational Field 


To address the problem of the motion of flux through in the halo, I consider a 
one dimensional diffusion problem of particles in a downward magnetic field 
of strength g with diffusion coefficient f3. 

The one dimensional Fokker-Planck equation in is 


dt 2 




(35) 


Let us start the problem at time t = 0 with f{z,0) = S{z). Dehne 
< {Azy/2 >) = f5. Then, with < Az >= —gt, the Fokker-Planck equation 
becomes. 


df{z,t) , ^df 


(36) 


Next transform this equation by setting z to z' = z + gt^{2 which gives 


df{z',t) 


dt 


dz'"^ 


whose solution is 


f{z\z) = 




■v/4^^ 

1 


exp 


: exp 


4:/3t 

{z + gty2)^ 


Now at any time t the number of particles with z > Ris 


(37) 


(38) 


m = 


f{z,t)dz 


'R 


\^Wt 
1 r 


exp 


(;^ + 9tV2y 


'R 


Af3t 


■dz 


du 


^ J{R+gt^/2)/y/{W) 


(39) 
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In other words, 


m 


Erie 


(R + gt‘l2)ly'Wt 


(40) 


where Erie is the compliment of the error function. 


This is Equation ([2]). 
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